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FOREWORD 


The  work  described  in  the  report  was  sponsored  by  the  Flight 
Dynamics  Branch  of  the  Air  Vehicle  Technology  Department  of  the 
Naval  Air  Development  Center.  Mr.  Robert  L.  Fortenbaugh,  Mr.  A. 
Piranian,  Mr.  Ronald  L.  Nave  and  Mr.  John  Clark  served  as  the  Navy 
Technical  Monitors  for  the  Naval  Air  Development  Center.  Mr.  Carmen 
J.  Mazza  was  NADC  program  manager. 

The  project  was  managed  by  Dr.  Harold  Stalford.  This  report 
is  a  summary  of  the  simulation  study  documented  in  Volume  II  and  a 
summary  of  the  real  flight  test  data  analysis  documented  in 
Volume  III. 
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1.  INTRODUCTION 

This  study  was  motivated  by  the  inability  of  empirical  and 
theoretical  techniques  to  accurately  predict  the  nonlinear  aerodynamics 
of  the  high  angle  of  attack  /sideslip  domain  for  aircraft.  The  Naval  Air 
Development  Center  (NADC)  took  the  direction  to  extract  the  nonlinear 
aerodynamics  from  actual  flight  test  data  using  system  identification 
and  modeling  techniques.  Success  in  this  direction  allows  predicted 
effects  to  be  correlated  with  the  actual  effects  to  determine  the  scaling, 
trending,  reformulation,  etc.  that  must  be  applied  to  the  prediction  to 
match  the  actual.  By  performing  this  correlation  for  several  existing 
aircraft  it  will  build  confidence  in  the  designer's  ability  to  predict  the 
performance  and  stability  and  control  characteristics  of  proposed  air¬ 
craft. 

The  T-2C  is  a  light  jet  trainer  aircraft  which  was  instrumented 
and  operated  by  NADC  for  the  express  purpose  of  producing  high  angle 
of  attack  flight  data.  The  data  along  with  a  wind  tunnel  model  of  the 
T-2C  was  supplied  to  Dynamics  Research  Corporation  (DRC)  for  pro¬ 
cessing  with  an  advanced  system  identification  method  called  the  Esti¬ 
mation-Before-Modeling  (EBM)  technique.  The  desired  results  of 
applying  the  EBM  technique  to  flight  data  are  (1)  optimal,  smoothed, 
time  history  estimates  of  the  vehicle  state  parameters  and  (2)  state - 
dependent  models  for  the  coefficients  of  the  forces  and  moments  acting 
on  the  vehicle. 

The  scope  of  this  aircraft  system  identification  study  is  confined 
1)  to  processing  16  maneuvers  of  synthetic  T-2C  flight  test  data,  which 
generated  using  the  T-2C  wind  tunnel  model  and  a  theoretical  predic¬ 
tion  model,  and  2)  to  processing  18  maneuvers  of  actual  T-2C  flight  test 
data. 
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1. 


The  objective  of  the  study  is  two-fold: 

Develop  and  refine  the  EBM  method  for  application  to  the 
identification  of  aircraft  aerodynamic  characteristics  in  the 
non-linear  (high  angle  of  attack/ sideslip)  flight  regime. 

2.  Apply  the  EBM  method  to  both  simulated  and  actual  flight 
test  data  for  a  T-2C  jet  trainer  aircraft  at  high  angles  of 
attack  and  sideslip. 

A  major  contribution  of  this  study  is  the  approach  taken  to  identify, 
from  actual  flight  test  data  in  high  angle  of  attack/  sideslip  flight  regimes, 
the  highly  nonlinear  aerodynamic  stability  and  control  characteristics. 

The  EBM  technique  is  a  two-step  approach  to  the  identification  of  aircraft 
aerodynamics  (See  Figure  1-1).  The  first  step  which  is  the  nonlinear 
estimation  part  uses  Gauss-Markov  processes  to  model  the  accelerations 
due  to  aerodynamic  forces,  (e.  g. ,  drag,  sideslip  and  lift)  and  moments 
(e.g. ,  roll,  pitch  and  yaw)  .  It  employs  an  extended  Kalman- Bucy  filter/ 
Bryson- Frazier  smoother  to  generate  smoothed  estimates  of  the  states, 
aerodynamic  forces  and  moments  time  histories,  and  measurement  biases 
and  scale  factors.  In  the  second  step  which  is  the  modeling  part,  the 
Stepwise  Multiple  Linear  Regression  (SMLR)  method  is  used  together  with 
the  subspace  modeling  to  identify  the  nonlinear  dependency  of  the  aerodynamic 
coefficients  on  aircraft  states  and  controls. 

Any  system  identification  method  must  perform  both  estimation 
and  modeling.  A  unique  feature  of  the  EBM  technique  is  the  separation 
of  the  estimation  of  the  aircraft  states  (i.  e. ,  the  first  step)  from  the 
modeling  of  the  aerodynamics  coefficients  (i.  e. ,  the  second  step).  It 
is  this  separation  that  gives  the  EBM  technique  its  power  as  a  modeling 
tool.  It  decouples  the  highly  complex  modeling  problem  into  six  simpler 
ones.  Instead  of  the  problem  of  model,  g  all  six  coefficients  simultaneously, 
each  coefficient  is  modeled  independently  of  the  others.  In  addition. 


OVERVIEW  OF  EBM  METHOD 


Aerodynamic 

Model 


instead  of  identifying  a  model  of  each  coefficient  for  each  maneuver, 
a  single  model  for  each  coefficient  is  identified  for  all  maneuvers. 
Furthermore,  the  modeling  of  a  coefficient  is  conducted  over  sub¬ 
spaces  of  the  state-control  domain.  The  models  over  the  subspaces 
are  low  order  models.  Synthesis  of  the  subspace  models  provides  the 
identified  nonlinear  model.  Thus,  a  highly  nonlinear  model  is  identified 
through  the  synthesis  of  identified  low  order  models. 

The  summary  of  the  application  of  the  EBM  technique  to  syn¬ 
thetic  T-2C  data  is  given  in  Section  2.  (Volume  II  contains  the  full  set 
of  results).  The  conditions  of  the  simulation  study  are  described  in 
Section  2. 1.  The  estimation  results  which  are  the  output  of  the  first 
step  of  the  EBM  technique  are  summarized  in  Section  2.2.  The  model¬ 
ing  results  which  are  the  output  of  the  second  step  of  the  EBM  technique 
are  described  in  Sections  2.3  and  2.4  for  lateral  and  longitudinal  co¬ 
efficients,  respectively.  Prediction  results  of  the  identified  model  are 
summarized  in  Section  2.  5. 

The  summary  of  the  application  of  the  EBM  technique  to  real 
T-2C  data  is  presented  in  Section  3.  (The  full  set  of  results  are  given 
in  Volume  III).  The  T-2C  flight  test  conditions  and  a  description  of  the 
maneuvers  are  summarized  in  Section  3.1.  A  summary  of  the  estima¬ 
tion  results  is  given  in  Section  3. 2;  these  results  are  the  output  of  the 
application  of  the  first  step  of  the  EBM  technique  to  the  real  T-2C  data. 

A  summary  of  the  identified  lateral  model  of  the  actual  T-2C  aircraft  is 
presented  in  Section  3.3.  The  longitudinal  modeling  results  are  sum¬ 
marized  in  Section  3.4. 

The  conclusions  and  recommendations  of  the  study  are  summar¬ 
ized  in  Section  4. 
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2.  SUMMARY  OF  SIMULATION  STUDY  RESULTS 
2. 1  CONDITIONS  OF  THE  SIMULATION  STUDY 

A  T-2C  wind  tunnel  data  model  and  a  theoretical  prediction  model 
of  the  dynamic  derivatives,  all  supplied  by  NADC,  formed  the  6  DOF 
model  used  in  the  simulation  study  to  generate  synthetic  test  flight 
data.  The  synthetic  maneuvers  were  made  as  realistic  as  possible  by 
using  the  control  input  time  histories  of  the  actual  T-2C  test  flight 
data.  Sixteen  maneuvers  of  available  T-2C  actual  flight  test  data  were 
used  to  provide  the  control  time  histories  and  the  initial  conditions  for 
generating  sixteen  synthetic  maneuvers.  These  are  summarized  in 
Table  2. 1-1.  The  designation  SF1M1  in  the  first  column  denotes  that 
the  controls  are  from  maneuver  1  of  flight  1  of  the  actual  flight  test 
data.  The  type  of  control  is  described  in  the  third  column  for  each 
maneuver.  Some  characteristics  of  the  synthetic  maneuvers  are  given 
in  the  last  column. 

Under  the  conditions  described  above,  the  equations  of  motion 
were  integrated  while  exciting  the  synthetic  T-2C  model  with  the  control 
inputs  from  the  actual  flight  test  data.  When  this  was  done,  it  was  found 
that  a  few  maneuvers  diverged  after  a  few  seconds,  and  the  responses 
for  most  maneuvers  were  totally  different  from  the  actual  flight  test  data. 
The  synthetic  model  exhibited  frequencies  and  amplitudes  that  were  very 
much  higher  than  the  actual  flight  test  data.  An  example  of  this  behavior 
is  shown  in  Figure  2. 1  of  Volume  II.  In  that  example,  the  actual  aircraft 
motion  had  a  maximum  angle  of  attack  of  25°,  but  the  synthetic  maneuver 
had  an  angle  of  attack  as  high  as  60°.  The  problem  was  found  to  be  in  the 
theoretical  prediction  model  of  Cn  .  By  trial  and  error,  the  constant 
value  of  -.06  for  Cn  was  found  t<?  give  stability  to  the  simulated  data  and 
to  provide  synthetic  Maneuvers  which  agreed  better  with  the  actual  flight 
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test  data.  Fixing  the  value  of  C  was  no  loss  to  checking  out  the 

"p 

EBM  technique  since  there  was  a  gain  in  stability  and  since  other 

dynamic  derivatives  C .  ,  C  „  ,  and  C  provided  similar  non- 

p  r  r 

linearities  as  exhibited  in  the  original  model  of  C  . 

np 

2.  2  ESTIMATION  RESULTS  -  FIRST  STEP  OF  THE  EBM  METHOD 

The  time  history  of  the  states  for  each  of  the  sixteen  synthetic 
maneuvers  were  fed  into  NADC  supplied  measurement  equations  to  pro¬ 
duce  noise-free  measurements.  The  measurements  were  corrupted 
with  white,  zero-mean  gaussian  noise  having  the  standard  deviations 
supplied  by  NADC.  These  standard  deviations  are  shown  in  Table  2.2-1 
under  the  heading  "Measurement  Noise  Level"  for  each  measured  vari¬ 
able.  The.  corrupted  measurements  were  processed  using  an  extended 
Kalman-Bucy  filter/Bryson-Frazier  smoother  to  obtain  estimated  values 
of  the  states  and  of  the  forces  and  moments.  For  the  purpose  of  checking 
the  performance  of  the  filter /smoother,  the  estimated  values  were  fed 
into  the  measurement  equations  and  the  RMS  values  of  the  difference  be¬ 
tween  these  "smoothed"  measurement  values  and  the  noise-free  meas¬ 
urements  were  computed.  The  RMS  values  averaged  for  all  sixteen 
maneuvers  are  presented  in  the  last  column  of  Table  2.2-1.  This  table 
shows  how  much  of  the  measurement  error  was  taken  out  by  the  estima¬ 
tion  process.  The  amount  of  noise  corrupting  the  synthetic  data  is  given 
in  the  fourth  column.  The  amount  of  noise  remaining  after  the  estimation 
process  is  shown  in  the  last  column.  For  example,  the  noise  level  on  the 
angle  of  attack  is  reduced  by  the  estimation  process  from  .  C028  rad  to 
.0005  rad,  which  is  a  reduction  of  82%  so  that  only  18%  of  the  original 
noise  remains  in  the  smoothed  data.  Noise  levels  are  reduced  by  at  least 
50%  for  most  of  the  other  measured  variables  as  well.  Using  Table  2.2-1 
we  make  the  following  observations  of  the  estimation  results.  The  estima¬ 
tion  process  has  reduced  the  error  level  from  .  24  to  .  12  for  the  accelera¬ 
tion  measurement  A  X  NOSE  in  the  x-direction.  The  errors  on  the 
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acceleration  measurements  in  the  z-direction  (i.e.,  AZCG,  . . .,  AZL) 
have  been  reduced  by  a  factor  of  3  or  4.  The  errors  on  the  acceleration 
measurement  in  the  y-direction  have  been  reduced  by  a  factor  of  3.  The 
errors  on  Q  and  0  have  been  reduced  by  a  factor  of  3.  The  errors  on 
the  roll,  pitch  and  yaw  rates  have  been  reduced  by  less  than  a  factor  of 
2.  The  errors  on  air  speed  are  less  by  a  factor  of  7.  The  errors  on 
a  and  0  have  been  cut  down  by  a  factor  of  5  or  6.  These  estimation 
results  demonstrate  that  the  Kalman-Bucy  filter /Bryson -Frazier 
smoother  performs  well  as  a  nonlinear  estimator  of  aircraft  states, 
forces  and  moments. 

2.  3  LONGITUDINAL  MODEL  IDENTIFICATION  RESULTS  - 

SECOND  STEP  OF  THE  EBM  TECHNIQUE 

The  smooth  data  provided  by  the  estimation  process  is  first 
analyzed  for  its  distribution  in  the  a,  0  state  space.  The  distribution 
of  the  smooth  data  for  the  range  [-15°,  -10°]  is  shown  in  Table  2.3-1. 
Each  number  in  a  subspace  is  the  number  of  times  that  the  sixteen  man¬ 
euvers  had  data  points  on  the  subspace.  About  fifteen  data  points  are 
needed  to  model  the  longitudinal  coefficients  on  a  subspace.  The  twenty- 
five  subspaces  containing  at  least  fifteen  points  are  enclosed  with  a  heavy 
black  line. 

Other  subdivisions  were  made  in  the  &e  space  to  define  the  final 
set  of  subspaces  for  modeling.  That  is,  the  longitudinal  coefficients 
were  modeled  on  subspaces  defined  by  ct,  and  6^.  The  subspaces 
are  described  in  Table  2.3-2.  The  number  of  subspaces  are  presented 
in  the  last  column  for  the  corresponding  and  0  intervals.  From 
the  table  we  see  that  the  a  ,  0  and  6  modeling  space  is  subdivided  into 
88  subspaces  to  be  modeled  with  the  SMLR  technique.  Instead  of  trying 
to  match  the  entire  data  with  one  equation  containing  a  very  large  number 
of  terms,  we  have  divided  the  data  from  the  16  maneuvers  into  88  groups 
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Table  2.  3-1  Density  of  Modeling  Data  for  the  Subspaces  Contained  in  the 

Region  Defined  by  |j3|*  0*  to  30  0  and  a  *  -4°  to  35°  and  by 

6  Tuvin  a  -15*  and  6  may  a  -10° 
e  e 


such  that  we  can  model  each  region  with  a  simple  low-order  model. 

Then  the  resulting  88  low-order  models  are  combined  into  a  global  non¬ 
linear  model. 

From  Table  2.3-2  we  observe  that  the  subspaces  containing  enough 
data  for  modeling  do  not  form  a  covering  of  the  entire  a,  /3,  4^  space  and 
there  are  holes  (e.  g.  third,  fourth  and  fifth  rows  of  the  table).  We  see  that 
the  controls  used  were  inadequate  to  generate  data  in  several  regions, 
especially  at  high  0  and  also  for  the  combination  of  low  a  and  high  4  . 
There  are  not  enough  data  points  for  a  beyond  35°. 


Identification  of  C  ,  C  and  C  using  the  SMLR  technique 
x  z  m 

was  performed  in  all  subspaces  and  then  synthesized  into  a  global 

model.  The  global  model  is  documented  in  Section  4.3  and  4.4  of 

Volume  II.  Models  of  C  ,  C  and  C  and  their  derivatives  with 

x  z  m 

respect  to  a,  0  and  6g  were  identified  for  0°sas35°,  0°s|(3|  s  10°, 

and  -24°  £  6  s  0°.  The  dynamic  derivatives  C  +  C  and  C  +  C 

e  z  z  -  m  m  • 

o  q  a  q  a 

were  identified  for  0  s  a  s  35  .  Examples  of  comparisons  between  the 

identified  models  and  the  synthetic  model  (i.  e. ,  wind  tunnel  and  theoretical 

prediction)  are  given  in  Figures  2.3-1,  2.3-2,  2. 3-3  and  2.  3-4.  For 

6e  =  -25°  and  |  $\  *  5°  the  identified  model  of  Cx  (o$,  0°sq!s35O,  is 

compared  with  the  wind  tunnel  model  in  Figure  2.3-1.  The  wind  tunnel 

values  are  represented  by  a  ’’square"  symbol  and  the  identified  values  are 

denoted  by  a  "triangular"  symbol.  The  identified  values  lie  almost  on  top 

of  the  wind  tunnel  values. 


oa.TR  b-U  DU  BETA  *  S  DEG 

t  a  UTUMTO) 


i  17.  SO  26.25  35 .00 

ALPHA! DEG) 


Figure  2.3-1  Comparison  between  the  identified  model  of  the  C 
coeffjcient  and  the  T-2C  wind  tunnel  model  for  X 
P  s  5  and  =  -25  :  Synthetic  data 
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For  5e  =  -15°  and  |/?|  =  10°  the  identified  model  of  C ^  (a), 

0°  s  a  s  35°,  is  compared  with  the  wind  tunnel  data  in  Figure  2.  3-2. 
All  identified  values  are  highly  accurate  except  for  a.  =  30°.  We  see 
from  Table  2.3-1  that  almost  no  data  points  were  available  for  modeling 
in  the  neighborhood  of  a  =  30°,  |0|  =  10°,  and  6e  =  -15°.  Yet,  the 
identified  value  at  a  -  30°  is  close  to  the  wind  tunnel  value. 


oq.tr  chi  ota 
90  *  UTIMTCO 


Figure  2.3-2  Comparison  between  the  identified  model  of  the  C 
coefficient  and  the  T-2C  wind  tunnel  data  for  2 
0  =  10  and  6  *  -15  :  Synthetic  data 

For  |  0|  ■  5°,  the  identified  model  of  C  (a),  0°sas  35°, 

o  ^  9 

6e  3  “25  ,  -15  ,  -10  ,  0  ,  is  compared  with  the  wind  tunnel  data  in 

Figure  2.3-3.  The  four  curves  represented  in  the  figure  correspond  to 
3  -25°,  -15°,  -10°  and  0°.  Observe  how  the  two  wind  tunnel  data 
curves  with  =  -15°  and  -25°  coincide  at  low  a  but  differ  consider¬ 
ably  at  high  a.  This  nonlinear  separation  in  the  wind  tunnel  model  was 
accurately  identified. 
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The  identified  model  of  C  (o)  +  C  (a),  0°  s  as  350  is 

m  m* 

q  a 

compared  in  Figure  2.3-4  with  the  theoretical  prediction  model  of 

Cm  +  .  We  have  once  again  excellent  agreement  between  the 

q  u 

identified  model  that  came  out  of  the  simulation  study  and  the  synthetic 
model  that  went  into  it. 


0  anmno 


Figure  2.3-4  Comparison  between  the  identified  model  of  the 

combined  derivatives  C  +  C  and  the  theoretical 

m  m* 

prediction  model:  Syntheti?  data  a 
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The  results  presented  in  Figures  2.3-1,  2.3-2,  2.3-3,  2.3-4 

are  typical  of  the  identified  models  for  C  ,  C  and  C  obtained  in 

x  z  m 

,  the  simulation  study. 

2.4  LATERAL  MODEL  IDENTIFICATION  RESULTS  -  SECOND 
STEP  OF  THE  EBM  TECHNIQUE 

The  Of  and  0  spaces  were  used  to  define  the  subspaces  for 
lateral  modeling.  The  Of-space  was  divided  into  18  intervals  and  the  0 
space  was  partitioned  into  3  intervals:  0°-2°,  2°-5°,  and  5°-10°. 

This  partitioning  on  a  and  0  produced  54  subspaces.  There  was 
enough  data  with  20°  s  a  s  25°  and  10°  s  |  0|  s  15°  to  provide  an 
additional  subspace,  making  a  total  of  55  subspaces  for  lateral  model¬ 
ing.  Thus,  instead  of  trying  to  find  one  very  complex  model  of  high 
order  for  the  lateral  coefficients,  we  reduced  the  modeling  problem  to 
finding  55  low -order  models  in  the  55  subspaces,  which  give  a  global 
nonlinear  model  when  synthesized. 

The  coefficients  C  ,  C .  and  C  were  identified  in  each  of  the 
y  *  n 

subspaces  as  low -order  functions  (i.  e.,  linear  or  quadratic  terms)  of 

the  variables  a,  0,  S  ,  6  ,  p,  r  and  0.  Synthesizing  the  identified 

models  produced  nonlinear  models  of  the  static,  control  and  dynamic 

derivatives.  In  particular,  for  the  range  of  a  and  0  with  0  s  30° 

and  Os  |j3|  s  10°  we  identified  nonlinear  models  for  the  following 

derivatives  of  the  coefficients  C  ,  C  „  and  C  : 

y  4  n 
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K 


? 


n 


c  (a,  0) 

YP 

C,  (a,  0) 

P 

C  (a,  0) 

P 

C  (a,  0) 

y6 

r 

C.  (a,  0) 

6 

r 

C  (a,  0) 

6 

r 

C  {a,  8) 

y6 

a 

C,  (a.  0) 

6 

a 

C_  (a,  (3) 

6 

a 

For  the  a- range  Osas  30°  we  identified  nonlinear  models  (except  for 

C  which  is  a 

n 

P 

coefficients  C 

y 

constant)  for  the  following  dynamic  derivatives  of  the 

,  C.  and  C  : 

*  n 

C 

-Z _ 

C 

n 

C  (a) 

Jr 

Cf  « 

r 

Cn  <*> 
r 

Ci  <“> 

P 

C  (c* 
n 

P 

(a,  0) 
a 

Cn  fa  P) 

u 

C£  (a,  0) 
rat 

Cn  (at.  0) 
ra 

C£  (0!.  0) 
pa 

The  identification  process  was  successful  in  separating  C  and  C  , 

nr  n0 

in  separating  and  C ^  ,  and  in  correctly  determining  that  the  deriva- 

r  0 

tives  C ,  ,  and  C  ,  are  insignificant.  (The  synthetic  model  values  for 

P  P 

these  two  derivatives  are  zero). 
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Some  typical  examples  of  the  identified  models  are  shown  in 
Figures  2.4-1,  2.4-2,  2. 4-3  and  2.4-4.  The  identified  model  of 
Cn  (o £,  0)  for  0  =  0°,  2^,  5°  and  10°  is  compared  with  the  wind  tunnel 
data  values  in  Figure  2.4-1.  For  0=2°  and  0  =  5°,  C n  initially 
decreases  in  magnitude  with  a,  increases  near  stall  (i.e.,  a  =  15°), 
decreases  rapidly  after  stall  becoming  nearly  zero  at  a  =  20°  and  then 
remains  very  small  in  magnitude.  For  0  =  10°,  decreases  in  mag¬ 
nitude  with  a,  falls  sharply  in  the  stall  region  and  dips  to  a  negative 

value  at  a  -  18°.  With  further  increase  in  a,  the  value  of  C  increases 

n 

slightly.  Very  good  agreement  is  indicated  between  the  identified  model  and 
the  wind  tunnel  model.  The  high  nonlinearities  present  around  the  stall 
region  were  accurately  identified. 

We  see  in  Figure  2.4-2  that  the  identification  process  obtained 

good  estimates  of  C  (a,  0)  at  0  =  2°.  It  is  positive  at  low  a,  declines 

n0 

in  value  rapidly  through  stall  and  becomes  extremely  small  or  negative  in 
the  post  stall  region. 

CM  S  a  ESTIMATED 


ALPHA 

Figure  2.4-2  Comparison  between  the  identified  model  of  the 
derivative  Cn  and  the  T-2C  wind  tunnel  data 

for  0  =  2°:  Simulation  Study 
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In  Figure  2.4-3  we  see  that  the  identified  values  of  the  dynamic 
derivative  C  is  nearly  a  constant  value  of  -.06  which  is  the  synthetic 

"p 

model  value  used  in  this  simulation  study. 
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Figure  2.4-3  Comparison  between  the  identified  model  of  the 
dynamic  derivative  Cn  and  the  synthetic  model 

of  C  :  Simulation  Stu ,Sy 


In  Figure  2.4-4  we  observe  the  comparison  between  the  identified 

model  for  C  A  and  the  theoretical  prediction  model  used  in  this  Simula - 
P 

tion  study.  The  roll  damping  derivative  C.  stays  constant  up  to  an 


I 


angle  of  attack  of  10°  when  it  starts  decreasing.  Immediately  after 
stall  it  jumps  to  an  unstable  value  of  0.3  and  stays  positive  up  to  an 
a  of  20°.  Beyond  that  the  roll  damping  is  very  small.  The  identified 
model  follows  this  highly  nonlinear  variation  extremely  well. 


a  OTOWTO 
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Figure  2.4-4  Comparison  between  the  identified  mod*il  of  the 

dynamic  derivative  C  ^  and  the  theoretical  prediction 

model:  Simulation  StudJ? 


2 . 5  PREDICTION  RESULTS 


The  identified  global  model  of  the  six  coefficients  C  ,  C  ,  C  , 

x  y  z 

C C  and  C  was  assessed  as  to  how  well  the  identified  model  pre- 
*  m  n 

diets  (1)  the  sixteen  maneuvers  used  in  the  identification  and  (2)  a 
maneuver  not  used  in  the  identification  process. 


For  the  first  case,  excellent  agreement  between  the  original 
synthetic  (i.e.,  "true")  response  and  the  response  predicted  by  the 
identified  model  was  obtained  for  most  of  the  maneuvers.  An  example 
of  the  comparison  of  the  responses  of  a  and  /3  for  maneuver  SF2M3 
is  given  in  Figure  2. 5-1.  The  excellent  matches  of  a  and  )3  are  typical 
for  the  other  states,  forces  and  moments. 


Figure  2.5-1  Comparison  between  the  predicted  and  the  true  time 
histories  of  a  and  0  for  SF2M3 


For  some  maneuvers  the  response  of  the  identified  model  went 
into  regions  where  no  models  were  identified.  Extrapolating  the  identi¬ 
fied  model  into  these  regions  led  to  a  response  that  was  sometimes  dif¬ 
ferent  from  the  "true"  response. 

For  the  second  case,  new  control  inputs  (Figure  2.5-2)  were 
chosen  to  excite  both  the  identified  model  and  the  synthetic  model. 
Sinusoidal  inputs  were  used.  The  responses  of  a.  and  0  are  compared 
in  Figure  2.5-3.  We  notice  that  the  response  predicted  by  the  identified 
model  and  the  actual  response  are  in  good  agreement. 


O 


Figure  2.5-2  The  controls  6  6  and  &  used  in  predicting  a  new 

maneuver 
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3.  SUMMARY  OF  ACTUAL  FLIGHT  TEST  DATA  RESULTS 


3. 1  FLIGHT  TEST  CONDITIONS  AND  DESCRIPTION  OF 
MANEUVERS 

The  T-2C  flight  test  data  were  gathered  with  gear  up  and  speed 
brake  retracted  under  the  following  conditions: 

True  Airspeed  -  75  to  150  meters /sec 

Mach  No.  -  0.25  to  0.4 

Altitude  -  6000  to  7600  meters 

The  throtiies  were  fixed  during  maneuvers.  The  Naval  Air  Development 
Center  provided  four  sets  of  flight  data  consisting  of  eighteen  maneuvers, 
amounting  to  over  600  seconds  of  20  Hertz  sampled  data.  The  type  of  input  used 
to  generate  each  maneuver  is  described  in  the  last  column  of  Table  3.1-1. 

The  first  column  gives  the  date  of  the  flight  set  and  the  second  column 
denotes  the  flight  number  and  the  maneuver  number  that  goes  with  the  in¬ 
put  used. 

The  maneuvers  F1M1,  F4M9  and  F6M1  have  360°  rolls  and  the 
maneuver  F4M10  has  two  360°  rolls.  The  angle  of  attack,  o i,  ranges 
from  -4°  to  40°  and  the  sideslip  angle,  /3,  covers  the  range  from  -14° 
to  26°.  Most  of  the  maneuvers  have  a  trim  a  between  10.5*  and  11.5° 
and  a  trim  V  (i.  e.,  true  airspeed)  between  90.0  and  95.0  meters/sec. 

The  maneuvers  F4M7  and  F4M8  have  a  trim  a  between  5.5°  and  6.5° 
and  a  trim  V  between  115  and  120  meters /sec. 

3.  2  ESTIMATION  RESULTS  -  FIRST  STEP  OF  THE  EBM  METHOD 

The  extended  Kalman-Bucy  filter /Bryson-Frazier  smoother  as 
discussed  in  Volume  H  was  used  to  process  the  eighteen  maneuvers 
F1M1,  F2M1, . . . ,  F6M1  to  obtain  the  estimation  results  of  this  study. 

Bias  states  and  scale  factor  states  were  added  to  the  filter  /smoother 
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to  handle  the  instrumentation  biases  and  scale  factors  that  are  inherent 
in  real  data.  The  estimates  of  the  scale  factors  and  the  biases  are 
given  in  Tables  3.2-1  and  3.2-2.  The  last  column  of  Table  3.2-2  lists 
the  flight  numbers  for  which  estimated  biases  were  observed.  Note 
that  no  biases  were  observed  for  oc  and  V  for  the  maneuvers  of  Flight 
Number  4.  We  see  from  Table  3. 1-1  that  the  maneuvers  of  Flight  No.  4 
were  run  on  August  26,  1976  and  that  the  other  maneuvers  were  run 
before  December  30,  1975.  The  roll  rate  and  pitch  rate  biases  have  large 
magnitudes,  and  there  was  no  consistency  in  the  estimates  of  the  sideslip 
angle  bias.  The  measured  variables  are  described  in  Table  3.2-3. 

The  estimated  values  of  the  states,  forces  and  moments  of  FI  Ml 
are  compared  with  their  corresponding  measured  values  in  Figures 
3.2-1,  3.  2-2  and  3.  2-3.  Comparisons  of  a  and  of  (3  are  given  in 
Figure  3.2-1.  Therein,  the  estimated  values  are  plotted  with  crosses 
at  one  second  intervals.  Note  how  the  estimator  smooths  out  the  measured 
seesaw  values  caused  by  the  vibration  of  the  noseboom.  In  Figure  3.  2-2 
observe  how  the  estimator  smooths  through  the  loss  of  the  airspeed 
measurement  shortly  before  8.0  seconds.  The  estimator  smoothed  out 
the  high  noise  level  on  the  yaw  rate.  The  estimation  results,  of  which 
those  of  F1M1  are  typical,  were  very  good  for  all  states,  forces  and 
moments  of  all  the  maneuvers  of  the  flights,  with  the  exception  of  F4M10. 

O 

In  F4M10  the  roll  angle  passes  through  180  several  times;  the  measure¬ 
ment  of  the  roll  angle  is  not  a  good  measurement  in  the  neighborhood  of 
180°. 
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Table  3.2-1  The  estimated  values  of  the  scale  factors  of  a  and  0 


Measured  Variable 

Estimated 

Scale  Factor 

a. 

.78 

8 

.805 

Table  3.2-2  The  estimated  values  of  the  biases  of  the  measured 
variables 


Measured  Variable 

Estimated 

Bias 

- 1 

Estimated  in  the 
Maneuvers  of  the 
Following  Flights 

0 £ 

O 

2 

1,  2,  3  and  6 

8 

Variable  up  to  5° 

an 

V 

-3.  9  m/sec 

1,  2,  3  and  6 

P 

-6. 1  /sec 

all 

q 

-2. 7°/sec 

an 

6 

a 

Fore-Aft  accel. 

0 

1.8 

,  2 

-0.  6  m/sec 

all 

an 

Normal  accel.  (c.g.) 

.  2 

-1. 6  m/sec 

aH 

Normal  accel.  (off c.g.) 

2 

1. 0  m/sec 

an 

Lateral  accel.  (c.g.) 

2 

-.  13  m/sec 

aU 

Rate  gyro  alignment 

O 

-3 

all 

Table  3.2-3  List  of  the  measured  variables  of  the  T-2C  actual 
flight  test  data 


! - 

Measured  Qty. 

Description  of  Measured  Quantity 

a  (alpha) 

angle  of  attack 

(3  (Beta) 

sideslip  angle 

V 

total  airspeed 

P 

roll  rate 

q 

pitch  rate 

r 

yaw  rate 

0  (Theta) 

pitch  angle 

<t>  (Phi) 

roll  angle 

z 

attitude 

AZ  CG 

vertical  acceleration  measured  at  c.g. 

AZR 

vertical  acceleration  measured  at  right  wing 

AZL 

vertical  acceleration  measured  at  left  wing 

AZ  Nose 

vertical  acceleration  measured  at  nose 

AZ  Tail 

vertical  acceleration  measured  at  tail 

AY  CG 

lateral  acceleration  measured  at  c.  g. 

AY  Nose 

lateral  acceleration  measured  at  nose 

AY  Tail 

lateral  acceleration  measured  at  tail 

AX  Nose 

Fore -aft  acceleration  measured  at  nose 

3.3  SUMMARY  OF  THE  LATERAL  MODEL  IDENTIFICATION 
RESULTS  -  SECOND  STEP  OF  THE  EBM  METHOD 


The  states  a  and  0  were  used  for  defining  the  subspaces 
for  modeling.  The  subspaces  in  a.  were  created  by  partitioning  a 

O  O 

at  one  degree  intervals  between  -4  and  25  and  at  five  degree  inter¬ 
vals  between  25°  and  40°.  The  subspaces  in  0  were  handled  by  using 
spline  functions  with  knots  at  five  degree  intervals. 

The  estimated  results  of  the  eighteen  maneuvers  F1M1, .... 

F6M1  were  processed  by  the  stepwise  multiple  linear  regression 

(SMLR)  technique,  the  modeling  phase  of  the  EBM  system  identification 

method.  Nonlinear  models  of  the  lateral  coefficients  C  ,  C,  and  C 

y  *  n 

were  identified  as  functions  of  the  variables- a,  0,  p,  r,  6  and  6  . 

a  r 

Table  3.3-1  contains  a  listing  of  the  identified  derivatives,  where  the 
limit  0max(od  is  the  maximum  value  achieved  by  the  sideslip  angle 
0  over  the  eighteen  maneuvers  at  the  angle  of  attack  a. 

Table  3.3-1  The  derivatives  of  the  lateral  coefficients  that  were 
identified  with  nonlinear  models,  -4  s  a  s  40  , 

0  s  |  0|  s  0max'(a) 


C  derivatives 

C.  derivatives 

C  derivatives 

y 

Xj 

n 

c  (4  0) 
y0 

C£  (a,  0) 

0 

C  (a,  0) 
n0 

- 

C,  (a) 

C  (a) 

nfi 

a 

a 

c  ta> 
yt 

C,  (a) 

C  (a) 

n6 

r 

r 

r 

- 

ci  <a> 

C  (a) 
n 

p 

P 

- 

Ct  (4) 

Cn  (4) 

r 

r 
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The  SMLR  technique  found  the  derivatives  with  respect  to  /3  to 

be  statistically  insignificant.  The  derivatives  C  (a)  and  C  (a)  were 

yp  yr 

masked  by  the  excessively  high  measurement  noise  level  on  the  yaw  rate 
measurement.  For  in  the  sideslip  equation  of  motion  the  yaw  rate  multi¬ 
plies  the  total  airspeed  whose  product  presents  a  high  noise  level  in  the 
identification  of  the  C  derivatives. 

y 

The  identified  state  and  control  derivatives  were  compared  with 
the  wind  tunnel  model  and  the  dynamic  derivatives  were  compared  with 
the  theoretical  prediction  models  of  Bihrle  Applied  Research.  Typical 
examples  of  these  comparisons  are  given  in  what  follows. 

The  comparison  of  the  C  (a,  0),  -4°  s  a  s  40°,  (3  =  9°  is 

yfl 

K  o 

shown  in  Figure  3.3-la.  The  curves  match  well  between  a  =  6  and 
a  *  35°.  Outside  of  this  a  range  where  there  was  sparse  data,  the 
identified  model  is  still  reasonably  good.  This  is  a  typical  example  of 


the  good  match  between  the  identified  model  of  C  and  the  wind  tunnel 
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This  same  derivative  is  also  shown  in  Figure  3. 3 -lb  for  the  range 
o  £  |0|  S  10°,  a  =  22.5°;  the  EBM  technique  identified  the  "deep  well". 


Figure  3.  3 -lb  Comparison  between  the  identified  model  and  the  wind 
tunnel  data  '  of  C  (a,  0)  at  a  =  22.  5°. 
y0 


The  identified  model  of  C.  (a)  is  given  in  Figure  3.3-2.  It 

\ 

a 

compares  well  with  the  wind  tunnel  data  at  a  below  stall,  and  its 
drop  in  controllability  above  stall  coincides  with  that  of  the  wind  tunnel 
data.  For  a  between  18  and  26,  the  identified  model  has  slightly 
smaller  magnitude  than  the  wind  tunnel,  but  they  coincide  above  an  a  of 


Figure  3.3-2  Comparison  between  the  identified  model  and  the  T-2C 
wind  tunnel  data  for  C. 

a 
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The  identified  model  of  the  dynamic  derivative  (a)  is 
presented  in  Figure  3.  3-3.  At  low  a  the  identified  mode? has 
smaller  values  than  the  theoretical  prediction  model.  The  two 
curves  cross  near  stall  and  the  identified  model,  hovers  about  the  zero 
axis  between  16°  and  20*  whereas  the  theoretical  prediction  goes  posi¬ 
tive  and  reaches  a  peak  of  0.30  per  rad.  The  two  curves  are  close 
for  a  above  20°. 


Figure  3.3-3  Comparison  between  the  identified  model  and 

the  theoretical  prediction  mBdel  of  Bihrle  Applied 
Research 

The  identified  model  of  the  dynamic  derivative  (o)  is 

r 

plotted  in  Figure  3.3-4.  It  has  a  slightly  steeper  slope  at  low  a 
and  drops  off  more  sharply  at  stall  than  the  theoretical  prediction 
model.  Above  a  =  18*  the  identified  model  lies  slightly  under  the 
theoretical  prediction. 


RLPHflf  OEQ) 


Comparison  between  the  identified  model  and 

the  theoretical  prediction  model  of  Bihrle  Applied 
Research 


The  identified  model  of  C 


(a,  (3;  for  (3  =  9  is  given  in 
closely  matches  the  wind  tunnel 
data  at  low  a,  where  it  is  nearly  a  constant  value.  At  stall  the 
identified  model  drops  off  with  the  wind  tunnel  value  where  they  take 
on  some  negative  values.  The  two  curves  have  close  agreement  above 


The  identified  ma 


flliPHfU.DEi 


Figure  3.3-5  Comparison  between  the  identified  model  and  the  T-2C 
wind  tunnel  data  for  C  for  (3  =  9° 


The  identified  dynamic  derivative  C  (a)  is  plotted  in  Figure 

np 

3.3-6  where  it  is  compared  with  the  theoretical  prediction  model  of 
Bihrle  Applied  Research.  The  identified  model  agrees  with  the  theoreti¬ 
cal  prediction  at  angles  of  attack  below  stall.  The  identified  model 
reaches  its  highest  positive  value  of  .  113  at  a  =  16. 5°  whereas  the  theo¬ 
retical  prediction  model  continues  to  grow  reaching  its  peak  of  0.  365  at 
a  =  18°.  After  reaching  a  peak  at  a=  16.5°  the  identified  model  drops 
off  to  a  value  of  -.076  at  oe  =  20.5°,  rises  to  .048  at  a  =  22.5;  xnd 
drops  to  -.051  at  a  =  23.5  before  leveling  off  about  -.075  over  the  a. 
range  [25°,  40°].  The  theoretical  prediction  model,  on  the  other  hand, 
drops  from  its  peak  at  a  =  24°,  then  rises  to  .  195  at  a  -  30°  before 
dropping  to  .06  at  ol  -  38°. 

Of  particular  interest  is  the  difference  between  these  two  models 
over  the  a  range  [20°,  35°].  In  this  range  the  theoretical  prediction 
model  has  values  greater  than  ,  10,  while  the  identified  model  has  values 


Figure  3.3-6  Comparison  between  the  C  identified  model  and  the 

n 

P 

theoretical  prediction  model  of  Bihrle  Applied  Research 
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around  -.075,  except  for  the  sharp  rise  and  decline  between  a  =  21° 

and  23°  .  This  difference  is  highly  significant  when  correlated  with 

the  findings  of  the  simulation  study.  Therein,  the  initial  conditions 

and  the  recorded  control  time  histories  of  the  actual  T-2C  flight 

data  were  used  to  drive  the  synthetic  T-2C  model  composed  of  the  wind 

tunnel  model  and  the  theoretical  prediction  model  of  Bihrle  Applied 

Research.  The  synthetic  responses  were  drastically  different  from 

the  actual  flight  data  until  a  change  in  the  theoretical  prediction 

model  of  Cn  was  made.  By  trial  and  error  a  constant  value  of  Cn 
P  P 

=  -.  06  was  found  to  give  stability  to  the  simulated  data.  Note  that  the 

identified  model  of  Cn  (about  -.  075)  obtained  by  processing  the 

P 

actual  T-2C  flight  data  agrees  with  the  modification  of  C  obtained 

np 

by  trial  and  error  in  the  simulation  study.  This  result  provides  strong 

evidence  that  the  Bihrle  model  for  C  is  inaccurate. 

P 

The  identified  model  of  Cn  (a)  is  plotted  in  Fugure  3.3-7. 

It  is  a  good  match  to  the  theoretical  prediction  model  of  Bihrle  Applied 
Research.  Around  stall  there  appears  to  be  a  difference  between  the 
curves  of  about  1°  shift  in  a  . 


Figure  3.3-7  Comparison  between  the  identified  model  and  the 
theoretical  prediction  model  of  Bihrle  Applied  Research 
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3.4  SUMMARY  OF  THE  LONGITUDINAL  MODEL  IDENTIFICATION 
RESULTS  -  SECOND  STEP  OF  THE  EBM  TECHNIQUE 

The  longitudinal  coefficients  C  ,  C  ,  and  C  are  functions 

x  z  m 

of  the  following  variables: 

a  -  angle  of  attack 

0  -  sideslip  angle 

q  -  pitch  rate 

5  -  elevator  deflection 

e 

• 

a  -  a  rate 

The  states  a,  0  and  6g  are  used  for  defining  the  subspaces  for  model¬ 
ing.  Subspace  modeling  is  conducted  in  two  steps  for  the  modeling  of 
Cm«  The  first  uses  a  partitioning  of  the  a  space  [-4°,  40°]  together 
with  spline  functions  to  handle  the  subspaces  of  the  )3  and  the  &e  space. 
The  subspaces  in  a  are  generated  by  partitioning  the  a  space  into  one 
degree  intervals  [a.,  a.  +  l°],  a  =  -4°,  -3°,  ...  ,  24°  and  five  degrees 
intervals  [o^,  +  5°],  <x  *  25°,  30°,  35°.  The  eighteen  estimation  data 

files  (i.  e.,  one  for  each  of  the  maneuvers  F1M1,  F2M1,  ....  F4M10 
and  F6M1)  generated  by  the  estimation  phase  of  the  EBM  system  identi¬ 
fication  method  are  transformed  into  thirty-two  new  estimation  data 
files  for  modeling  by  using  the  partitions  on  a  and  by  collecting  together 
the  data  from  all  eighteen  maneuvers  having  an  angle  of  attack  within  a 
given  a -interval.  In  this  first  step  the  subspaces  in  j3  and  are 
handled  by  using  spline  functions.  Since  a  longitudinal  coefficient  is 
symmetric  with  respect  to  the  sideslip  angle,  a  quadratic  spline  func¬ 
tion  as  a  function  of  |  j3|  is  used. 

The  second  step  uses  a  partitioning  of  the  6  space  [-28°,  2°] 
together  with  spline  functions  to  handle  the  subspaces  of  a.  The  space 
of  is  partitioned  into  one  degree  intervals  ,  6  +  l°],  * 

-28°,  -27°,  ...  ,  1°.  The  eighteen  estimation  data1, files  generated  by 
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the  estimation  phase  of  the  EBM  method  are  transformed  into  thirty 
new  estimation  data  files  for  modeling  by  using  the  partitions  on 
and  by  collecting  together  the  data  from  all  eighteen  maneuvers  having 
an  elevator  control  deflection  within  a  given  6e-interval.  In  this  sec¬ 
ond  step  the  subspaces  in  a  are  handled  by  using  spline  functions. 

Nonlinear  models  of  the  C  coefficient  were  identified  for  the 

m 

following  derivatives: 

C  (a,  6  )  ,  -28°  s  6  s  1°  ,  a  MIN  (6  )  s  a  s  a  MAX  (6  ) 

me*  e  e  e 

a 

C  (a,  0)  ,  -4°  s  a  s  40°  ,  Os  |0|  s  0  MAX  (a) 

mu 

C  (a)  ,  -4°  s  a  £  40° 
m 

q 

C  (a)  ,  -4°  s  a  s  40° 

m  * 
a 

where  a  MIN  (6e>  is  the  minimum  value  of  the  angle  of  attack  for  all 
maneuvers  at  the  elevator  setting  and,  similarly,  a  MAX  (6^)  is 
the  maximum  value. 

The  identified  dynamic  derivative  C  (a)  is  compared  with 

mq 

the  theoretical  prediction  model  of  Bihrle  Applied  Research  in  Figure 

3.4-1.  The  identified  model  of  C  varies  between  -11.4  and  -8.0 

m 

q 

at  a  below  4.5°.  It  then  climbs  slowly  to  its  peak  of  5.61  at  a  =  13.5° 
before  falling  to  -7.  7  at  a  =  18.5°  where  it  levels  off  to  values  between 
-7.7  and  -6.3  for  a  beyond  18.5°.  The  theoretical  model  holds  almost 
level  with  a  value  of  -10.  3  below  a  =  14°.  It  climbs  sharply  to  a  peak 
of  6.  9  at  a  *  1.  6#  before  dropping  slowly  to  -4. 5  at  a  =  30°.  The 
identified  model  has  positive  values  for  a  between  11°  and  16°,  whereas 
the  theoretical  model  has  positive  values  for  a  between  15.  5  0  and 
22.5°  . 
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Figure  3.4-1  Comparison  between  the  C  identified  model  and 


the  theoretical  prediction  of  Bihrle  Applied  Research 
Actual  Flight  Test  Data 


An  example  of  the  identified  model  of  C  (a)  at  6  =  - 

ma  e 

and  its  comparison  with  the  wind  tunnel  is  gi/en  in  Figure  3.4-2 


Figure  3.4-2  Comparison  between  the  Cm  identified  model  and  the 

wind  tunnel  data  for  5 g  =  -*55.5°:  Actual  Flight  Test 
Data 
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The  estimated  values  of  C^(t)  were  compared  with  the  wind 

tunnel  values.  The  elevator  control  <>  was  used  to  subdivide  the 

e 

Cz(t)  estimates  obtained  from  the  eighteen  maneuvers  into  data  sub¬ 
sets.  These  subsets  were  ordered  with  respect  to  alpha  and  then 
plotted.  The  figures  are  contained  in  Volume  III  and  they  show  a 
comparison  with  the  wind  tunnel  data.  The  data  in  these  subsets 

were  never  processed  to  generate  an  equation  model  of  C  ,  but  the 

z 

figures  show  that  the  model  of  the  actual  aircraft  is  close  to  that  given 


by  the  wind  tunnel  analysis.  In  the  figures,  the  real  data  values  below 
stall  almost  coincide  with  the  wind  tunnel  values  which  implies  that  the 


actual  C  (a)  is  close  to  the  wind  tunnel  data.  Around  stall  there 


appears  to  be  a  hysteresis  effect  in  the  real  data  that  is  not  accounted 
for  by  the  wind  tunnel  model.  For  example,  this  can  be  observed  in 
Figure  3.4-3  for  6  =  -14.5°. 


Figure  3.4-3  Comparison  between  the  real  data  values  and  the  wind 
tunnel  data  for  C  at  6  =  -14.5° 


=  -14.5 


The  estimated  values  of  Cx  (t)  were  compared  with  the  wind 

tunnel  values.  In  figures  in  Volume  III  the  aerodynamic  coefficient  Cx 

is  shown  as  a  function  of  alpha  for  fixed  one  degree  intervals  of  the 

elevator  control.  The  range  of  the  elevator  control  covers  -28°  to  0°. 

The  figures  were  obtained  as  follows:  A  one-degree  6^-interval  was 

selected.  All  eighteen  maneuvers  were  searched  to  find  all  Cx  (t) 

estimates  having  a  corresponding  5  within  the  selected  4  -interval. 

e  e 

The  Cx  (t)  estimates  were  then  ordered  according  to  alpha  and  plotted 
as  a  "real  data  values"  curve  in  Volume  III.  The  estimated  states  were 
used  to  evaluate  the  wind  tun  el  data  in  order  to  get  a  comparison. 
These  curves  have  a  dependency  on  /3,  q  and  a.  built  into  them  by  the 
nature  of  the  maneuvers.  The  curves  exhibit  some  differences  between 
the  wind  tunnel  data  and  the  actual  aircraft.  A  hysteresis  effect  can 
be  observed  in  the  figures.  The  real  data  values  are  above  the  wind 
tunnel  values  at  low  alpha  and  at  high  alpha.  Near  stall  they  are  below 
the  wind  tunnel  values.  For  example,  this  can  be  observed  in  Figure 
3.4-4. 


4.  CONCLUSIONS  AND  RECOMMENDATIONS 
4. 1  CONCLUSIONS  OF  THE  SIMULATION  STUDY 

The  EBM  methodology  has  been  validated  under  controlled  but 
realistic  simulated  conditions.  The  T-2C  wind  tunnel  data  and  a 
theoretical  prediction  model  of  the  dynamic  derivatives  were  used  to 
generate  synthetic  data.  The  controls,  initial  conditions,  masses  and 
thrusts  of  sixteen  actual  T_2C  flight  test  data  records  were  used  to 
excite  the  T-2C  wind  tunnel  data  model  to  generate  the  sixteen  synthetic 
maneuvers.  The  synthetic  data  were  corrupted  with  realistic  measure¬ 
ment  noise  levels  and  then  processed  with  an  extended  Kalman- Bucy/ 
Bryson-Frazier  smoother  to  produce  excellent  estimated  values  of  the 
states  and  of  the  forces  and  moments.  Subspace  modeling  together  with 
a  step-wise  multiple  linear  regression  technique  were  used  to  identify  a 
global  state /control  dependent  model  of  the  force  and  moment  coefficients. 

The  identified  model  agrees  well  with  the  T-2C  wind  tunnel  data.  The 
high  nonlinearities  were  accurately  identified.  The  identified  global  model 
gave  excellent  predicted  responses  for  a  new  maneuver  not  used  in  the 
identification  process. 

4.  2  CONCLUSIONS  OF  THE  REAL  DATA  ANALYSIS 

The  EBM  methodology  was  successfully  tested  on  flight  test  data 
and  provided  good  aerodynamic  models  of  the  T-2C  aircraft. 

Within  the  first  step  of  the  EBM  method,  the  data  from  eighteen  maneu- 
ers  which  had  been  collected  over  a  period  of  ten  months  were  processed, 
and  large  biases  were  accurately  estimated  in  angle  of  attack,  sideslip 
angle,  total  airspeed,  roll  rate,  pitch  rate,  fore-aft  acceleration,  normal 
acceleration  and  rate  gyro  alignment.  The  method  satisfactorily  handled  the 
high  noise  levels  on  the  yaw  rate,  the  rudder  deflection,  the  total  airspeed, 
the  acceleration  measurements  (as  well  as  the  losses  of  some  measurements 
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over  extended  intervals  of  flight),  and  the  nose  boom  vibration  present 
in  the  angle  of  attack  and  sideslip  angle  measurements.  The 
estimation  results  which  are  documented  in  Volume  III  and  exemplified 
in  Section  3.  2  demonstrate  that  the  first  step  of  the  EBM  method  per¬ 
formed  well  on  the  T-2C  flight  test  data. 


Within  the  second  step  of  the  EBM  method  detailed  detailed 

nonlinear  models  of  the  aerodynamic  derivatives  were  identified  and 

compared  with  the  wind  tunnel  data.  .  Most  of  the  identified  state  and 

control  derivatives  matched  well  with  the  wind  tunnel  data.  The 

identified  dynamic  derivatives  were  compared  with  the  theoretical 

prediction  models  of  Bihrle  Applied  Research.  The  identified 

models  of  C  .  and  agreed  fairly  well  with  the  theoretical 
n  I* 

prediction  models  with  the  exception  that  the  identified  model  hovers 

about  a  zero  value  around  stall  while  the  theoretical  model  shows 

positive  values.  The  identified  model  of  C ^  drops  off  more  sharply 

r 

after  stall  than  the  theoretical  model.  The  identified  C  compares  well 

“r 

for  all  alpha,  but  C  only  matches  at  low  alpha.  At  high  alpha,  C 
is  predicted  to  have  Values  above  .10  but  the  identified  model  showed 
that  the  values  are  around  -.075.,  The  simulation  study  documented 
in  Volume  II  found  that  the  theoretical  prediction  of  Cn  is  inaccurate 
at  high  ot  and  that  a  constant  value  of  -.06  provided  a  Setter  model  for 


generating  synthetic  responses  which  mimic  more  closely  the  real 

data.  The  identified  model  of  Cn  obtained  by  processing  the  actual 

T-2C  flight  data  agrees  with  the  modification  of  the  theoretical 

model  of  C  obtained  by  trial  and  error  in  the  simulation  study. 
nP 

The  modeling  results  of  the  real  data  analysis  study  demonstrate 
that  the  EBM  system  identification  method  provides  accurate  nonlinear 
modeling  of  high  aJ  0  aerodynamic  stability  and  control  characteristics 
from  flight  test  data. 
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4.3  MAJOR  CONCLUSIONS  OF  THE  STUDY 

The  major  conclusions  of  this  study  are  summarized  as: 

•  Validated  System  Identification  Tool 

The  EBM  methodology  has  been  validated  in  great 
detail  on  both  synthetic  and  actual  flight  test  data 
and  has  been  shown  to  produce  accurate  nonlinear 
models.  The  EBM  technique  has  been  shown  to 
be  a  practical  tool  for  the  modeling  and  analysis 
of  vehicles  aerodynamics  that  works  on  test 
data  gathered  under  actual  flight  conditions. 

•  Advantage  of  Subspace  Modeling 

/ 

Subspace  modeling  is  an  inherent  feature  of  the  EBM 
technique  which  circumvents  the  problem  of  identi¬ 
fying  in  one  step  complex  nonlinear  models  over  the 
entire  «,  ^parameter  region  covered  by  the  flight 
test.  Instead,  the  problem  is  reduced  to  an  easier 
procedure  of  identifying  lower  order  models  within 
the  subspaces  and  then  synthesizing  a  global  model 
from  these  subspace  models.  In  addition,  breaking 
the  modeling  into  subspaces  provides  a  means  of 
determining  how  much  data  is  available  in  each  sub- 
space,  and  this  in  turn,  provides  a  measure  of 
model  accuracy  and  confidence  and  provides  an 
approach  for  determing  when  further  flight  testing 
is  required  to  fill  sparse  subspaces. 

•  Data  Versus  Modeling  Problems 

The  two  steps  of  the  EBM  technique,  estimation  in 
the  time  domain  followed  by  modeling  in  the  state 
domain,  were  established  specifically  to  enhance  the 
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process  of  distinguishing  data  problems  from 
modeling  problems.  Data  problems,  if  they  exist, 
are  uncovered  in  the  estimation  phase.  That  this 
works  is  exemplified  by  the  results  obtained  by 
the  application  of  the  EBM  technique  to  the  actual 
T-2C  data:  measurement  and  axis  alignment 
biases,  scale  factors,  etc. 


4. 4  RECOMMENDATIONS 

We  recommend  the  following  for  improving  future  flight  test 
data  for  system  identification  studies. 

•  I.  Instrumentation 

A.  Measure  the  Euler  yaw  angle  0.  This  will  permit  better 
accuracy  in  estimating  the  bias  on  the  yaw  rate  r. 

B.  Use  sensors  that  measure  the  absolute  values  of  the  Euler 
pitch  and  bank  angles,  0  and  0,  rather  than  measurements 
that  are  relative  to  arbitrary  Initial  settings.  This  will  allow 
the  use  of  all  maneuvers  to  estimate  the  biases  of  9  and  0, 
and  it  will  provide  a  check  on  the  condition  9  *a  at  trim. 

C.  Measure  critical  engine  parameters  for  estimating  thrust 
of  engines.  This  will  permit  using  the  actual  values  of  a  vary¬ 
ing  thrust  to  identify  a  better  drag  model. 

D.  Use  bank  angle  sensors  that  measure  accurately  through 
rolls  greater  than  180°.  This  will  permit  better  estimation  re¬ 
sults  at  these  bank  angles  during  high  roll  rates  and  it  will 
provide  continuity  through  0  ■  180°. 

E.  Measure  reliably  the  angular  acceleration  rates  p,  q  and 
r.  This  will  provide  a  check  on  off-c.g.  acceleration  measure¬ 
ments  and  more  accurate  estimation  of  C  (t),  C  (t)  and  C  (t). 
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F.  Use  nose  boom  with  less  vibrational  characteristics  than  that 
on  the  T-2C.  This  will  provide  better  a  and  0  data  for  estima¬ 
tion  and  modeling. 

G.  Use  airspeed  sensors  which  do  not  give  a  loss  of  measure¬ 
ments  over  extended  intervals  of  flight.  This  permits  better 
estimation  and  modeling  results. 

H.  Use  single  channels  for  all  measured  variables.  This 
provides  the  best  signal -to -noise  ratio  for  a  measured  quantity 
as  compared  to  using  a  commutator. 

I.  Use  pulse  coded  modulation  (PCM)  rather  than  pulse  ampli¬ 
tude  modulation  (PAM).  This  will  provide  cleaner  signals. 

Rigorous  Calibration  and  Consistency  Checkout 

A.  Calibrate  instrumentation  using  static  data  gathered  in 
hangar  (i.  e. ,  with  aircraft  in  hangar,  isolated  from  wind, 
turn  instrumentation  power  on  and  record  100  seconds  of 
data).  This  will  permit  the  isolation  and  removal  of  unneces¬ 
sary  biases  in  linear  and  angular  accelerations,  in  p,  q  and 
r,  and  in  6  and  0. 

B.  Calibrate  instrumentation  using  dynamic  data  gathered  on  a  long 
runway  (i.  e. ,  record  data  with  aircraft  jetting  down  runway  in 
takeoff  phase).  This  will  permit  the  isolation  and  removel  of 
unnecessary  biases  in  a  and  0. 

C.  Calibrate  instrumentation  using  dynamic  data  gathered  in 
trim  flight  for  various  trim  values  of  sideslip  angle  and  total 
airspeed.  This  will  permit  the  isolation  and  removal  of  biases 
in  the  controls  for  equilibrium  flight. 


D.  Check  out  the  consistency  of  the  measured  variables  »«ing 
several  representative  maneuvers  of  flight  data.  Use  the  equa¬ 
tions  of  motion  to  test  for  consistency.  This  will  provide  a 
checkout  on  scale  factors  and  axis  alignment  as  well  as  the 
compatibility  of  the  measured  variables. 

E.  Fly  calibration  maneuvers  before  and  after  each  flight 
set  of  maneuvers. 

HI.  Improve  Design  of  Maneuvers  for  System  Identification 

A.  Fly  more  maneuvers  and  vary  the  variables  such  as 
trim  airspeed  and  sideslip  angle. 

B.  Design  the  maneuvers  so  that  a  good  amount  of  data 
is  provided  at  high  and  low  a  and  at  high  0  extremes 
rather  than  a  few  isolated  data  points. 

C.  Start  and  end  each  maneuver  with  long  trim  times. 

IV.  Precise  Calculation  of  Aircraft  Constants 

A.  Compute  the  aircraft  mass  for  the  beginning  and  the  end 
of  each  maneuver. 

B.  Compute  the  limits  on  the  changes  in  the  moments  of 
inertia  due  to  weight  distribution  changes  such  as  fuel 
consumption  during  maneuvers. 

C.  Provide  calculations  on  the  location  of  the  center  of 
gravity  with  respect  to  the  aerodynamic  center  for  each 
maneuver. 
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